A New Monomeric Interpretation of Intrinsic 
Optical Bistability Observed in Yb 3+ - Doped 

Bromide Materials 



F. Ciccarello a 

A. Napoli a 
A. Messina a 
S.R. Luthi b 

a INFM and MIUR, Dipartimento di Scienze Fisiche ed Astronomiche 
dell'Universita di Palermo, via Archirafi 36, 90123 Palermo, Italy 

h Gemfire Corporation, 1220 Page Avenue, Fremont, CA - 94538, USA 



Abstract 

We present a mechanism able to show intrinsic bistable behaviour involving single 
Yb 3+ ions embedded into bromide lattices, in which intrinsic optical bistability 
(IOB) has been observed. The mechanism is based on the experimentally found 
coupling between the Yb 3+ ion and the totally symmetric local mode of vibration 
of the [YbBro\ 3 ~ coordination unit. The model reproduces the IOB observed in 
CsCdBr% : 1 % Yb 3+ and allows to understand the experimentally found presence 
of the phenomenon in the other bromides, but its absence in C s 3 Lu2Cl§ : Yb 3+ . 
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A physical system without mirrors of an external cavity, exhibiting two stable 
outputs for a single input intensity of an exciting laser, is said to exhibit in- 
trinsic optical bistability (IOB) [1,2,3]. Hehlen et al. [4,5,6,7] have discovered 
that Yb 3+ - doped Cs 3 X 2 Br 9 [X = Y, Lu] and CsCdBr 3 show IOB when 
exciting into the 2 F 7/2 (0) ^ 2 F 5/2 (0') and 2 F 7/2 (0) ^ 2 F 5 / 2 (2') electronic 
transitions of the Yb 3+ ion with a narrow band laser of suitable near-infrared 
(NIR) frequency. On the contrary, no IOB was observed under excitation into 
the 2 -F"7/ 2 (0) — > 2 -F5/ 2 (l') transition [8,9]. When doped in these crystalline bro- 
mide materials some Yb 3+ ions form strongly interacting Yb 3+ — Yb 3+ dimers 



Email address: messina@fisica.unipa.it (A. Messina). 



Preprint submitted to Elsevier Science 



2 February 2008 



[10,11]. Due to this interaction, excitation of the Yb 3+ 2 F 7 / 2 — > 2 -F5/2 transi- 
tion (around 10000 cm" 1 ) results in both, 2 F 5 / 2 —* 2 F 7 / 2 emission of single 
Yb 3+ ions in the NIR spectral range, as well as Yb 3+ — Yb 3+ dimer emission 
in the visible (VIS) at twice the resonance frequency of the single Yb 3+ ion 
[6,12,13] through a cooperative luminescence mechanism [14]. The IOB phe- 
nomenon consists of a temperature and excitation power dependent hysteretic 
behaviour of the transmitted, NIR and upconverted VIS emission intensity 
(with the system switching output intensity at two critical incident excitation 
powers). To date, it is discussed if the phenomenon is due to the presence of 
strongly interacting Yb 3+ — Yb 3+ dimers or if it is a single Yb 3+ ion prop- 
erty. The theories relying on the former interpretation are usually referred to 
as "dimeric" while those supporting the latter idea as "monomeric". Dimeric 
theories substantially rely on the existence of intra-pair electronic coupling of 
Yb 3+ — Yb 3+ dimers [4,5,15,16]. According to these theories, IOB should be 
observed a fortiori in Yb 3+ - doped Cs 3 Lu 2 Cl 9 , as this material is isostruc- 
tural to Yb 3+ - doped Cs^X 2 Br 9 [X = Y, Lu] [17,18] and exhibits a stronger 
electronic exchange interaction for Yb 3+ — Yb 3+ dimers [10]. However, no ex- 
perimental evidence for IOB was found in V6 3+ -doped Cs 3 Lu 2 Cl 9 [8,9]. Fur- 
thermore, in Cs^X 2 Br 9 : Yb 3+ (X = Y, Lu) only a few Yb 3+ ions form dimers 
at low concentration of impurity ions. Nevertheless, both bromide lattices ex- 
hibit IOB even at a low Yb 3+ - doping level [4,5,6]. A monomeric mechanism 
to account for IOB in these systems has been proposed by Gamelin et al. [9] 
who presented a phenomenological model, according to which IOB depends 
on a single - Yb 3+ ion property in these materials. In this model, IOB is at- 
tributed to the interplay of a laser heating effect and a strongly increasing, 
non-linear dependence of the material's absorbance with the internal temper- 
ature. With its phenomenological nature, this theory succeeds in reproducing 
many experimental trends, but does not provide the microscopic mechanism 
of IOB. In this letter, we present the first microscopic, monomeric mechanism 
that allows a single active ion to exhibit all essential aspects observed in IOB. 
To date, the microscopic theories presented in literature describing how a col- 
lection of two-level atoms can exhibit IOB are based on electronic coupling 
between two or more near absorbers [3,5,15,19]. In these cases, the electronic 
interaction between the coupled absorbers provides the feedback loop required 
to obtain intrinsic bistability. Together with the intrinsic two-level atom non 
linearity and the use of single-atom density matrix approach (an indispens- 
able step in order to obtain non linear equations), this feedback loop gives 
rise to inversion dependent atomic resonance frequency renormalization, and 
thus possible IOB. With a few exceptions, the role played by the host lattice, 
and in particular of the first coordination sphere of the Yb 3+ , in the mecha- 
nism responsible for IOB in V6 3+ -doped bromide host materials has generally 
been regarded as secondary. However, Hehlen et al. measured different inter- 
nal temperatures for the system on the two branches of the hysteresis loop 
[5]. This suggests that with the first switch occurring, the lattice undergoes 
heating. Additionally, Liithi [8] has found that Cs 3 Lu 2 Br g : 10% Yb 3+ shows 
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bistable behaviour even when excited with a chopped laser and off-times as 
long as 75 msec. This means that the system prepared on the upper branch 
of the bistability region does not decay to the lower branch for times as long 
as 75 msec after turning off the exciting laser beam. If IOB were due to a 
cooperative mechanism based on electronic Yb 3+ — Yb 3+ coupling, and the 
feedback loop provided by the same two ytterbium ions, it would be hard to 
justify this behavior. After 75 msec, all ions in their excited state are expected 
to have decayed to the ground state, and there seems to be no physical reason 
for the system to save memory of its past history. The fact that the nuclear re- 
laxation times of the lattice are much longer than the electronic ones, suggests 
that the lattice nuclear degrees of freedom play a relevant role in the bista- 
bility mechanism. Supported by the findings of a Raman spectroscopy study 
on CsCdBr 3 : 1% Yb 3+ [20] showing strong coupling between the electronic 
2 F 7 / 2 (0) F 5 / 2 (2') transition of Yb 3+ and the totally symmetric vibrational 
local mode of the [YbBr 6 ] 3 ~ coordination unit (into which Yb 3+ is embedded 
in all bromide lattices showing IOB), we propose that this vibrational cou- 
pling provides the indispensable feedback loop to observe IOB (involving only 
a single Yb 3+ ion). 

Our physical system consists of the electronic structure of a single Yb 3+ ion 
and the totally symmetric vibrational local mode A\ g of the \YbBre] 3 ~ octa- 
hedron. We treat the former as a two- level system (\a) and \b) ground and 
excited states, respectively, hu energy difference, \i) (j\ = i,j = a,b) and 
the latter as a quantum harmonic oscillator (a and a + annihilation and cre- 
ation operators, respectively, and angular frequency ujk). We write the free 
Hamiltonian H of the system as 

H = H Y b + H k = huj cr bb + hu k a + a (1) 

For the interaction Hamiltonian Vyb-k between the ion and the A\ g local 
mode we choose a bilinear operatorial form (in order for the oscillation centre 
to depend on the electronic state in which the ion is [21]) 

V Yb -k = € (a bb - a aa ) (a + a + ) (2) 

When a coherent classical electromagnetic field is applied to the ion, using 
the rotating wave and the electric dipole approximations, the coupling be- 
tween the Yb 3+ ion and the external field is described by the observable 
V Yb -L = ^-((Tbae~ iuJt + cr o6 e ia '*), where Q and u stand for the Rabi and the 
laser frequency. Denoting the total Hamiltonian operator of the system with 
H = H + Vyb-k + Vyb-L, and taking into account both the ion spontaneous 
emission and the dissipative interaction of the ion and the local mode with 
the bulk modes of the lattice (described by a suitable superoperator S), the 
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dynamics of the system obeys the generalized Von-Neumann equation 

F = - t -[H,F] + SF (3) 



where F is the density operator of the whole system. F(t) can be written as 
F(t) = Tr k F{t) Tr Yb F(t) + F Yb - k (t), where Tr k F and Tr Yb F are the reduced 
density operators of the ion and the local mode, respectively, while F Yb ~ k {t) 
represents the correlations existing at the time t between Yb and k. Assuming 
Vyb-k <C H , we approximate F(t) with its factorized part 

F(t) ~ Tr k F(t) Tr Yb F(t) (4) 



This is the main hypothesis of our model. In this way, analogously to other 
theories [3,5,19,22], we assume that single-ion density matrix approach is al- 
lowed. A necessary condition for this approach to be valid is that Vyb-k is much 
smaller than the free Hamiltonian H (later on we will show this is the case). 
It must be remarked that the feasibility of this treatment - indispensable to 
have non linear equations and thus possible IOB - poses a delicate and lively 
debated [15,23,24] theoretical problem questioning even the existence of a mi- 
croscopic explanation of some IOB phenomena [23]. We justify assumption 
(4) with the ability of our model to reproduce the observed IOB phenomena. 
Inserting eq.(4) in eq.(3) and doing the trace over k and Yb, we obtain two 
coupled master equations in Tr k F(t) and TrYbF {t). Performing the variable 
change p' ab = p ab e luJt , it can be proved that the ion dynamic equations in p bb 
and p' ab differ from the usual optical Bloch equations in the rotating frame for 
the presence of an effective resonance frequency ufj^ (t) = uoq + 2|(a + a + ){t), 
linearly dependent on the mean value of the local mode operator (a + a + ). 
This proves that the dynamic response of the ion to the applied field depends 
on the state of the local mode (in its turn depending on the ion state). It is easy 
to show that for stationary conditions (when setting all temporal derivatives 
to zero) the resonance frequency can be written as 

U)l fS = UJq- (3(p bb - p aa ) (5) 



where 



with 7 being the relaxation rate of the local mode. We thus find that the overall 
effect of the coupling between the electronic transition of the Yb 3+ ion and the 
Ai g local mode of the [YbBr & } 3 ~ unit, consists of an ion resonance frequency 
renormalization linearly dependent on inversion (p bb — p aa ). Introducing the 
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absorption parameter s = 2p b b (0 < s < 1), it is possible to prove that, under 
stationary conditions, s and the incident excitation power P must obey the 
equation 

c 3 s 3 + c 2 s 2 + ( Cl + bP) s -bP = (7) 



with 




r, = I J — ) C 2 = -8 — -p— Ci 

3e 2 F ab 2 




cos y? — — (8) 



nmce hu ^iTh 

where is the homogeneous linewidth of the electronic transition, 8 = uj — uo 
is the detuning of the laser frequency, F ab = (2mu /3h)\(a\r\b)\ 2 stands for 
the oscillator strength of the transition \a) — > |6), ip is the angle between 
the laser polarization versor and (a|r|6), and Ti is the ion longitudinal decay 
time. Eq.(7) is a third degree equation in the unknown s dependent on many 
parameters. It can be shown that the condition for the existence of a P-range 
inside which, for any value of P, two stable solutions for s are found (IOB) is 

l>"4 (9) 
1 h z 



This provides a threshold value of the local mode - Yb 3+ ion interaction 
strength required to observe intrinsic bistability compatibly with our micro- 
scopic feedback mechanism. Moreover, eq.(7) shows that, for an increasing 
value of the bistability parameter /3/Th, the incident power bistability range 
moves toward higher values of P and, in addition, becomes wider. We now use 
our theory to reproduce the two incident powers at which switching occurs 
measured in CsCdBr 3 : 1 % Yb 3+ [7]. Guillot-Noel et al. successfully per- 
formed an analogous simulation to test their dimeric theory [15] . In our study 
we use: 5 = (perfect resonance), u = 10602.8 cm' 1 [7], T h = 0.6 cm' 1 [7], 
1/Ti = 7.8 • 10~ 4 s [13], F ab = 10~ 6 and, additionally, we take cos 2 ip = 1/2, 
assuming that the incident laser beam is unpolarized [5]. Using the bistability 
parameter to optimise our simulation, good agreement is achieved between 
the theoretical prediction and the experimentally measured values [7] for the 
two incident powers at which switching occurs, as shown in Fig. 1. Using the 
optimised value of f3/T h = 7.5, the ion-local mode coupling constant e [eq.(2)] 
can be estimated. From eq.(6) and taking into account the fact that the lat- 
tice nuclear relaxation times are typically very long compared to l/ou^, we 
get e ~ \\/f3 oo k . With oo k = 191 cm' 1 [20], we obtain e = 14.66 cm' 1 for 
CsCdBr 3 : 1% Yb 3+ . This value confirms our a priori assumption that Vyb-k 
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Fig. 1. Absorption parameter s of the Yb 3+ ion as a function of the incident intensity 
P (kW/cm 2 ) in CsCdBr 3 : 1% Yb 3+ at 7 K calculated from eq. (7) (fl/T h = 7.5, 
S = 0, uj = 10602.8 cm' 1 , T h = 0.6 cm" 1 , 1/Ti = 7.8 • W~ A sec, F ab = 10~ 6 ). 
The dashed line represents unstable solutions. The grey shaded area represents the 
experimentally observed bistable region for this compound [7]. 



is much weaker than Ho. With this value for e, the displacement of the oscil- 
lation centre for the local mode k, for an Yb 3+ ion undergoing the transition 
from | a) to |6), as implied by the operatorial form assumed for Vyb-k [ ec L-(2)], 
turns out to be in the order of 1 % . This result agrees with typical values 
found for lattice point defects due to incorporation of rare earth impurity ions 
[21], and confirms that the coupling constant e we get from our model is of 
the correct order of magnitude. 

Another relevant test for our theory is its ability to predict the presence of 
IOB in Cs 3 X 2 Br 9 : Yb 3+ [X = Y, Lit], but the absence of the phenomenon in 
Cs 3 Lu 2 Cl 9 : Yb 3+ [8,9]. In the bromide compounds, the angular frequency u k 
is 190 cm" 1 [12], similar to the one found in Cs 3 Er 2 Br 9 [25]. From comparison 
of the absorption and emission spectra, we estimate Th and e to be comparable 
to the ones measured in CsCdBr 3 : Yb 3+ [8,20]. These assumptions made, we 
expect the bistability parameter (3/Y h for Cs 3 X 2 Br 9 : Yb 3+ [X = Y, Lu) to 
have a comparable value to the one calculated for CsCdBr 3 : 1% Yb 3+ , and 
can understand why IOB is found in the former compounds. For the chloride, 
ujk is 285 cm -1 [25]. From inspection of the absorption and emission spectra, 
we estimate to fall in the range 0.3 — 0.7 cmT x [8], and the Yb 3+ - A lg 
local mode coupling constant e to be about 5 times smaller than in bromide 



6 



compounds. Based on these values, the bistability parameter /3 /Th is estimated 
to be lower than the critical threshold of 3y / 3/2 required for IOB. We finally 
notice that absorption spectra of Cs 3 Lu2Brg : Yb 3+ show that the Yb 3+ ion 
2 F 7 / 2 (0) F 5 / 2 (l') transition exhibits a very weak coupling to the A lg local 
mode compared to — 0' and — 2'. This explains why no IOB is observed 
under excitation into this peculiar transition. 

In summary, we have presented the first microscopic mechanism according 
to which a single Yb 3+ ion is able to exhibit IOB when doped in Cs 3 X 2 Br$ 
[X = Y, Lu] or CsCdBr 3 . The mechanism couples the electronic two- level 
system of the Yb 3+ ion with the totally symmetric local mode of the vibration 
of the [YbBro\ 3 ~ coordination unit. The theory has been used to reproduce 
basic aspects of the IOB observed in CsCdBr 3 : 1% Yb 3+ . We emphasize 
that our theory not only allows to understand, why IOB has been observed 
in CsCdBr 3 : 1% Yb 3+ and Cs 3 X 2 Br 9 [X = Y, Lu], but also, why it is 
absent in the chloride compound Cs 3 Lu2Clg : Yb 3+ . These results provide 
new insight into the microscopic mechanism underlying the IOB phenomenon 
and guidelines for expanding the material basis in the search for IOB systems, 
and their potential applications. 

We gratefully acknowledge inspiring discussions with O. Guillot-Noel of EN- 
SCP, Paris and M.P. Hehlen of M. Lujan Jr. of Neutron Scattering Centre, 
Los Alamos National Laboratory, Los Alamos, New Mexico. 

Note: Octahedral coordination of the [YbBro\ 3 ~ unit, instead of the actual C 3 
symmetry, is assumed throughout the text with no impact on our conclusions. 
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